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Abstract: New arguments are presented to emphasize the interest of the infrared finite 
couphng approach to power corrections in the context of Sudakov resummation. The more 
regular infrared behavior of some pecuhar combinations of Sudakov anomalous dimensions, 
free of Landau singularities at large A^/, is pointed out. A general conflict between the in- 
frared finite coupling and infrared renormalon approaches to power corrections is explained, 
and a possible resolution is proposed, which makes use of the arbitrariness of the choice of 
exponentiated constant terms. A simple ansatz for a 'universal' non-perturbative Sudakov 
effective coupling at large Nf follows naturally from these considerations. In this last ver- 
sion, a new result is presented: the striking emergence of an infrared finite perturbative 
effective coupling in the Drell-Yan process at large Nf (at odds with the infrared renor- 
malon argument) within the framework of Sudakov resummation for eikonal cross sections 
of Laenen, Sterman and Vogelsang. Some suggestions for phenomenology at finite Nf, 
alternative to the shape function approach, are given. 
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1. Introduction 

Tlie notion of an infrared (IR) finite coupling, and the related concept of universality, to 
parametrize power corrections in QCD has attracted much attention for a long time Q, 
Q. In the present note, which is an extended version of a talk given at the FRIF workshop 
on non-perturbative effects in jets (Paris, January 10-14 2006), I display further evidence 
in favor of this assumption in the more specific framework of Sudakov resummation. A 
more detailed account of the topics covered here is postponed to a future publication. The 
paper is organized as follows. Section 2 gives a phenomenological incentive for an IR finite 
effective coupling to parametrize the tail of the Sudakov peak for the thrust distribution. 
Section 3 gives a theoretical incentive, where some remarkable IR properties of specific 
combinations of Sudakov anomalous dimensions are pointed out. Section 4 discusses a 
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general conflict between the renormalon and the IR finite coupling approaches to power 
corrections, and indicates a possible resolution based on the arbitrariness of the choice of 
constant terms exponentiated in Sudakov resummation along with the large logarithms. 
This issue is further discussed in sections 5-7, and a closed form solution is given at large 
Nf in the case of deep inelastic scattering (DIS) in section 7.1. An application to the 
Drell-Yan process is given in section 7.2, where it is shown that the exponentiation of 
0{N^) terms advocated in naturally leads to the emergence of the simplest ansatz for 
an IR finite perturbative coupling (albeit being at odds with the IR renormalon argument). 
Alternatively, an attempt to reconcile the IR renormalon and IR coupling approaches leads 
to another ansatz for a 'universal', but non-perturbative, Sudakov effective coupling at 
large A'^^ proposed in section 8. Section 9 deals further with the issue of resummation 
of constant terms, and sketches a procedure for phenomenology at finite Nf. Section 10 
contains the conclusions. 



2. Incentive for IR finite coupling approach 

Let us take as an example the case of thrust, and define 



1 f^™^ , , . da 

/ or exp(— i/r) — = exp \E[Q 

c^tot Jo dr 



(2.1) 



where r = 1 — T. Consider now the case where v is large, which corresponds to r small, 
and assume the standard M exponentiation formula 



E{Q^, u) = S{Q^, v) + F(Q2) + o(xiy) 



(2.2) 



where H{Q'^) is a power series in as{Q'^) = as((5^)/47r with coefficients independent of v, 
and the 'Sudakov exponent' S{Q'^,v) is given by 



/■I r l-xQ'^ jl2 

5(Q2, u)= — [exp(-z.x) - 1] 4Cf / ^ Asik^) - 3CfBs{xQ 

Jo X J^2Q2 



(2.3) 



where ^^(A;^) and Bs{k^), the 'Sudakov effective couplings', should be considered as two 
physical 'effective charges' ^ which can be expanded as power series in ag^k"^) 



As{k') = as{k') + Aiaiik') + A2ai{k') + 



(2.4) 



and similarly for Bg(k'^). It is convenient to interchange the x and A;^ integrations, and 
write eg. (p. 3D as a sum of two 'renormalon integrals' B 



5(Q2, u) = ACf ^FAik'/Q^ u)As{k') - ?,Cf ^Fsik'/Q^ u)Bs{k''), (2.5) 
where F^ik'^ / , i^) and Fsik'^ /Q'^,^'), the 'Sudakov distribution functions', are given by 
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fk/Q ^™ 

FAik'/Q^ v)= — [exp(-i.x) - 1] , (2.6) 



and 



I note that 



FB{k'/Q\ v) = eM-^k'/Q') - 1. (2.7) 



FA{k'/Q',i^)^ln{k/Q) (2.8) 



and 



FBik^Q^i^)^-! (2.9) 

for z/ — > oo. Assume now that the effective coupUngs As{k'^) and Bs{k^) are IR finite, 
and setting As{k'^) = Air + AAsik^) and Bsik"^) ^ Bir + ABs{P), assume further that 
AAsiP) and ABs{k^) are 0(^2) foj. ^ q_ rj.^^^^ taking the — > oo hmit under the 
integrals, one gets 

S{Q\iy) - -ACFKi{u)AjR-3CFKo{iy)BiR 

+ ACfJ ^Hk/Q) AAs{e) + 3CF j ^ABsik^), (2.10) 

where the integrals over AAs{k'^) and ABs{k'^) are z^-independent (and IR convergent), 
whereas 

/■I 

Ki{iy) = — [exp(-z/x) - 1] Inx ~ 1/2 In^ (2.11) 
Jo X 

and 

Ko{iy) = / — [exp(-i/2;) - 1] ~ -Inz^ (2.12) 

JO X 

for 1/ — > oo. Thus we obtain at large u 

E{Q\ u) ~ -ACfKi{v)Air - ?,CfKo{u)Bir + C{Q^), (2.13) 

where 



C(Q2) = ^ ln(fc/Q) AAsik'') + 3C^ ^ABs{k'') + H{Q'') (2.14) 

is z/-independent. Hence for i/ — > oo 



/ exp(-i/r) — ^ U{Q^) exp r-2CFA7Rln2 i/ + 0(lnz/)l , (2.15) 
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where AA(Q^) is a non-perturbative Q-dependent normalization constant (its Q-dependence 
can eventually be estimated at large Q). Upon taking the inverse Laplace transform the 
present ansatz can thus potentially (if Aj^ > 0) reproduce the expected tail of the thrust 
distribution for r ^ 0, which is parametrized essentially by the IR value Ajr of the effective 
coupling As{k'^) {Bm gives a subdominant contribution) which could be fitted from the 
data. I note that no assumption that Q is large has been made, so this prediction goes 
beyond short distance physics. To fit the whole thrust distribution at all values of r, one 
should make a more complete ansatz for the IR behavior of As{k'^) and Bs{k'^)- This 
step requires a priori no more free parameters than the 'shape function' approach and 
provides an alternative to it. 



3. IR behavior of the Sudakov effective couphng in the Nf ^ oo Hmit 

Let us take as a case study the example of the scaling violation for the non-singlet structure 
function in deep inelastic scattering (DIS). Prom the standard exponentiation formulas in 
Mellin N- space (I shall adopt in general the notations of |9|), one gets immediately the 
relation at large for the 'physical anomalous dimension' § din F2{Q'^ , N)/dlnQ'^ 



d^^F^iQ'^^) _^CfH{Q^)+ACf f'dz '"'' ^ As[{l-z)Q']+0{l/N), (3.1) 

Jo i- — z 



where H(Q^) is again given as a power series in as(Q^) with iV-independent coefficients 
and 

AC^Asik') = A{a,{k^)) + (3.2) 

A (the universal 'cusp' anomalous dimension) and B are the standard Sudakov anomalous 
dimensions relevant to DIS, given as power series in a^: A{as) = A^ag + Aio?g + .. (with 
Aq = ^Cp) and B{as) = BqQs + Bia^ + ... I assumed the conjecture of relative to the 
vanishing to all orders of the third standard anomalous dimension D (the 'soft function'), 
is correct. As = + O(a^) shall be refered to as the 'Sudakov effective coupling', and we 
shall see (in the large Nf limit) that it has drastically different infrared properties then the 
individual Sudakov anomalous dimensions it is composed of. 

Large Nf analysis: the Borel transform i?[^5](ti) of the Sudakov effective coupling has 
been computed at large Nj. It is defined by 

As{k'^) = j exp j^-uln^J S[^5](u), (3.3) 



and is was found [11, M (this result is checked below using a different method) that at 



large Nf, with the 'nai've non-abelization' recipie |12] 



^ r . 1 / N 1 / , N sin vrii 1 / 1 1 , , , , 

B[As]{u) = -eM-du) - + , 3.4 

PO TTU 2, \\ — U \ — Ull' 
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where /?o = ^C2(G) - |A^/ is the one-loop coefficient of the beta function and d is a 
scheme-dependent constant related to the renormalization of fermion loops: d = —5/3 in 
the MS scheme and d = in the so called 'V-scheme' (or 'single dressed gluon scheme'). 
In the following I shall use the 'V-scheme' for simplicity. One should note that B[Ag]{u) 
is free of renormalons singularities, and in fact the corresponding standard perturbative 
series of ^4^ in powers of has a finite convergence radius. 

One first observes that for k = A (the Landau pole of the one-loop V-scheme coupling) 
the Borel integral in eq. (|3.3| ) converges and Ag{k'^ = A"^) is well-behaved and finite. This 
behavior is in striking contrast with that of the cusp anomalous dimension, which has also 



been computed ||15|| at large Nf 



A{as) = ^Cf 



sinvr/Joas r(4 + 2/3oas 



(3.5) 



TTf3o 6r(2 + /3oa,)2' 

where PqUs = l/ln/c^/A^ at large Nj. It is easy to see, since as oo for A; ^ A, that the 
ratio of gamma functions in eq.( |3.5D blows up in this limit, while the sin factor generates 
wild oscillations, resulting in a completely unphysical behavior around k = A. The much 
more tamed behavior of As{k^) suggests that it is well-behaved in the infrared region. It 
is actually possible to get an analytic expression for As{k'^), valid at all k"^. One finds 



with 



A 



simple / 7 2 
1 



1 , . Ei{t + m) — Ei{t — m) 
(A:^) + — ^exp(-t)^ 



27r/?o 



exp(-2t) 



Ei{2t + 2i7r) - Ei{2t - 2m) 



2i 



11 / , 
arctan(t/7r) 

2 vr 



(3.6) 



(3.7) 



where t = In/c^/A^ and Ei{x) is the exponential integral function. Eq.( |3.6| ) can equivalenly 
be written in term of the incomplete gamma function r(0, x) = ^ exp(— z) as 



Asik^) = 4*™^'^(A;2) + _i_exp(-t) 



27r/?o 



r(0, -t + iir) -r(0,-t -ivr) 



2i 



1 



27r/3o 



exp(-2t) 



r(0, -2t + 2m) - r(0, -2t - 2m) 
2i 



+ 7r 



To investigate the infrared behavior, I now make use of the fact that r(0, x) 
1/x + ...) for |x| oo, which implies 



exp(— x) 

X 



(3.8) 
(1- 



exp(— t) 



exp(-2t) 



r(o, -t + m) - r(o, -t - m)' 
2i 

"r(0, -2t + 2m) - r(0, -2t - 2m) ' 
2i 



vr 

t2 



27r 
'4i2 



(3.9) 
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for |t| — > oo. In the ultraviolet (UV) region {t — > +00) one thus obtain, as expected, a 
power series in as{k'^) = 1/Pot 

As{k') ^ asiP) + ^/3o alik') + - y)/?^ al{k') + .... (3.10) 

I note that Ag"^^'''^{k'^) behaves as + 0{a^) in this limit, while the last two terms on the 
right hand side of eq. (|3.8D are O(a^). A^^"^^''^ [k"^) turns out to coincide for real space-like 
k^ = fi"^ > with the (integrated) time- like discontinuity at k^ = —fi^ < of the one-loop 
coupling. I stress however that here it should be considered as an analytic contribution 
to the large Nj Sudakov effective coupling, an Euclidean coupling defined in the whole 
complex k"^ plane. The analytic continuation to complex k^ is conveniently provided by 
the identity: 

1 1 , , . 1 ln(t + in) - ln(t - in) 

arctanfr/vrj = . (-J-ll) 

2 TT TT 2i 

On the other hand, in the IR region (t — s- —00) A'^^"^^'''^ {k"^) reaches a finite limit of l//3o- 
Since the two incomplete gamma function contributions are 0{l/t'^), and thus also vanish 
for t — > —00, the full As{k'^) would have reached the same IR fixed point value of l//3oi 
were it not for the two IR divergent contributions 7rexp(— t) and 7rexp(— 2t). Thus in fact 
As{k'^) approaches an infinite (and negative) IR fixed point for fc^ —>■ 0: 

. .,2^ 1 1 A2 , , 

This is not by itself an unphysical behavior, except for the negative sign in the infrared, 
which cannot reproduce a vanishing Sudakov tail. Moreover As{k'^) is a causal function 
and has no unphysical Landau singularities on the whole first sheet of the k'^ plane: the 
only possible branch points occur for t = ±in, i.e. at k"^ / h? = — 1 on the time-like 
axis. The trouble is in the too strongly divergent IR behavior, which gives a divergent 
contribution ^0( (^-z)^ ) z — > 1 to the integral on the right-hand side of eq.(tj.l|). This 
finding makes it more plausible the speculation Q that there exists a non-perturbative 
modification 5As^Np{k'^) of the coupling at very small momenta which might turn the too 
strong IR fixed point of perturbative origin into a genuinely non-perturbative, but softer 
(eventually finite) IR fixed point, yielding an IR convergent Sudakov integral. A simple 
ansatz for dAg^j^pik"^)^ localized in the IR region, and which thus does not induce large 
OPE- violating |l| ultraviolet 0(1/(5^) power corrections, is given by 



1 


\ A2 / 






' k^\ 


Wo 











(3.13) 



Then 



As^Npik") = Asik") + 5As,Np{k' 



(3.14) 



^The straightforward modification of the coupling obtained by just subtracting the IR divergent contri- 
bution eq.(3.12) yields, when viewed from the UV side, a too large short distance 0{1/Q^) correction. 
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which indeed yields for k'^ a finite and positive IR fixed point: A^^Npik'^) — l//3o + 

4. A clash between the infrared finite couphng and the IR renormalons 
approaches to power corrections 

There is a potential clash between the infrared finite coupling and the IR renormalons 
approaches to power corrections (very closely connected to the well-known issue |jl^, |l^ 
of 1/Q corrections in Drell-Yan) which can be summarized in general terms as follows. 
Consider a typical 'renormalon integral' 

S{Q\N) = —F{k^/Q\N)As{k^), (4.1) 
and introduce its ('RS invariant' |jl6[) Borel representation by 

S{Q'^,N) = J^ du exp(^-uln^^ B[S]{u,N). (4.2) 
Using eq.(|3.3D one obtains the Borel transform of S as 



B[S]iu,N) = B[As]{u) F{u,N), (4.3) 

with 

/■I dx 

F(u,N)= — F(x,N) exp(-ulnx). (4.4) 
Jo X 

The factorized form | p^ , p^ ] of the expression should be noted. Suppose now the Borel 
transform B[As]{u) of the effective coupling Ag has a zero at some position u = uq, 
which is not shared by B[S]{u, N). Then necessarily F{u,N) in eq.( |4.3| ) must have a 
pole^ at u = Uq. If Uq > 0, this means an IR renormalon, and F{x,N) contains an 0(x"°) 
contribution for x ^ 0. Since this renormalon is not present in B[S]{u, N), the standard IR 
renormalon philosophy would conclude that no corresponding power correction is present. 
On the other hand, this power correction is still expected in the IR finite coupling approach, 
where As{k'^) is assumed to have a finite IR fixed point and the low energy part of the 
integral in eq.( [4.lD is well-defined. Then the distribution function F{k'^/Q'^,N) can be 
expanded in powers of A;^ for k"^ 0, yielding a non-vanishing power correction for each 
term in its low energy expansion, parametrized by a low energy moment of the effective 
coupling As{k^). Even if B[S]{u,N) is singular, rather than non-vanishing and finite, at 
u = Uq, there is still a clash, since the singularity of F{u,N) is necessarily stronger than 
that of B[S]{u,N) in presence of a zero of B[As]{u); thus the IR finite coupling approach 
will predict a more enhanced power correction than indicated by the renormalon argument. 



^This argument is usually presented |14j as the existence of a zero in the Borel transform Bfyls] causing 
the vanishing of the residue of a would-be renormalon contained in F{u, TV). Thus this remark is a particular 
case of the observation in that the renormalon residue is an all-order quantity, requiring the knowlegde 
of all the perturbative coefficients of the effective coupling As- 
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At this stage, there is no way to decide which of these two philosophies is correct without 
further information. I end this section by giving two examples of this situation at large 

i) DIS: there eq.(3^) implies that vanishes for any integer u, u > 3. 

ii) Drell-Yan: in this case the Sudakov effective coupling which occurs in the 
derivative of the Drell-Yan cross section is given by^ 

8CpAs,DY{k') = 2Aia,ik')) + (4.5) 



and the corresponding Borel transform at Nf = oo is |14| (dropping an inessential exp(cu) 
factor which corresponds simply to a change of scale): 

1 1 

B[As,DY]{n) = (4.6) 

/3or(l + ti)r(l/2-n) 

which vanishes for positive half-integer u. 

In both cases, the zeroes are absent from the corresponding Sudakov exponent, and the 
IR finite coupling approach will lead to the prediction of extra power corrections in the 
exponent {N^/Q^, N'^/Q^, etc.. .in DIS, and N/Q, N^/Q^, etc.. .in Drell-Yan). 



5. Ambiguities in Sudakov resummation and exponentiation of constant 
terms 

To make progress, we have to understand better the origin of zeroes in the Borel transform 
-B[A5](u) of the Sudakov effective coupling. In this section, I shall show it is possible to 
make compatible the renormalon and IR finite coupling approaches. The main point which 
will be developped is the following: there is an ambiguity in the choice of the Sudakov 
distribution function (or, alternatively, in that of the Sudakov effective coupling), which 
corresponds to the freedom to select an arbitrary set of 'constant terms' to be included in 
the Sudakov exponent. One can eventually make use of this freedom to reconcile the two 
conflicting approaches. 

To see this, let us consider the specific case of DIS scattering eq. (|3.1|) . Putting 

1-1 N-l _ 1 

S{Q\N)= dz— As[{l-z)Q% (5.1) 

JO I — z 

it is again convenient to write eq.(p.l|) as a renormalon integral (eq.( |4.1| )), where the Su- 
dakov distribution function is now given by 

F(A;VQ', N) = {1- k^/QY'^ - 1- (5.2) 

I first note that order by order the perturbative series of S{Q'^,N) contain both 0{N^) 
'constant terms' and terms which vanish as — > oo, e.g. 

^The natural occurence of this combination of anomalous dimensions has been pointed out independently 
from a different point of view by G.Sterman ||l7|]. 
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5(Q2, N) = (701L + 7oo)as(Q') + (7i2i>' + luL + 7io)as(Q') + - + 0{1/N), (5.3) 



where L = ln.N. 

A very useful simplification is achieved by making use of the following important 
scaling property: at large A^, F{k'^/Q'^,N) becomes a function of the variable e = 
only^, reflecting the relevance fl^ of the 'soft scale' Q'^/N. Indeed defining 



Fik^Q^, N) = G(e, iV) = (1 - e/N f-^ - 1, (5.4) 
and taking the N ^ 00 limit with e fixed one gets a finite result 

G{e,N) ^ G{e, 00) = G{e), (5.5) 

with 

G(e) = exp(-e) - 1. (5.6) 

Let us now redefine the Sudakov exponent by using G(e) as the new distribution function, 
i.e. eq. ([4.1[) is replaced by 

rQ'^ .11,2 

Sstan{Q\N) = —G{Nk^/Q^)As{k'), (5.7) 

where the index ^stan' stands for 'standard'. This step is legitimate since, order by order 
in perturbation theory, SgtaniQ'^ , Af) and S{Q'^,N) differ only by terms which vanish as 
N —I- 00, and thus share the same InN and constant terms (eq.( |5.3| )). Indeed, eq.(|5.7D can 
be written identically as 

Sstan{Q\N)= f^^ /Xp[-iV(l-z)]-l ^ , 

Jo 



and the equivalence between eq. (|5.1| ) and eq. (|5.8D up to 0{1/N) terms was proved in |19| . 
I now assume an ansatz of the form of eq. (|5.7D 



SneUQ^N) = J^ —GneUNkyQ^)A^^^{k^), (5.9) 

and show that a unique solution for Gnewi^k'^ /Q"^) and A^^'^^k'^) exists, under the condi- 
tion to reproduce all divergent InA^ terms, together with an (a priori arbitrary) given set 
of constant terms, i.e. 



SnevjiQ^N) = {j,,L + ^i;^^)as{Q') + {7i2L^ + 7uL + j'^rWs{Q') + ■■. + 0{1/N). (5.10) 

This statement can be checked order by order in perturbation theory. As is well-known, the 
IniV structure in eq.(p.3|) or ( [5. 101 ) arises from performing the integrals in eq.(p.7|) or ( ^.S| ) 

*For Drell-Yan, the corresponding scaling variable is eoY ~ ^y- 
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over the renormalization group ln(A;^/(5^) logs which show up when As{k'^) is expanded in 
powers of as{Q'^): 



Asik^) = asiQ^) + i-l3oHkyQ^) + Ai)a^,{Q^) + ... (5.11 

(where one should set fc^ = (1 — z)Q'^ if one uses instead eq. (|5.lD or ( |5.8D ). Using eq.(^ 
one thus gets 



SstaniQ^ N) = Ko{N)as{Q^) + [-poKiiN) + AiKo{N)] a^.iQ^) + (5.12) 

with 

^p(^) = J -^G{Nk^/Q^)\nP{k^/Q^) 

= / — G(e)lnP(e/A^). (5.13) 
Jo e 

Thus for large A'' one gets (see also eq.( |2.12| ) and ( p.llj )) 

Ko{N)= / -G{e) = coiL + coo + 0{l/N), (5.14) 
Jo ^ 



and 

rN 



de 

Ki{N)= / —G(e)ln(e/A^) = 012^2 + C11L + C10 + 0(l/iV), (5.15) 
Jo f 

where the Cj/s are known. Let us now try to modify G{e) — > G{e) + AG(e) = Gnewi^), in 
such a way the coefficients of all positive powers of IniV remain unchanged. I show that 
this cannot be achieved without changing simultaneously the Sudakov effective charge Ag. 
At 0{as) one gets cqo cqo + Acqo = Cqq^ , and one should require that 

r°° de 

Acoo = / — AG(e) < 00 (5.16) 
Jo ^ 

in order to preserve the leading log coefficient cqi. However at 0{af) one gets cn — > 
cii + Acii = with 

r°° de 

Acii = - / -AG(e) = -Acoo. (5.17) 







Thus the single log term in Ki{N) must he changed, and according to eq.( ^.12 ), this implies 



a correlated change Ai ^ Ai + AAi = Ai^^ if one wants to preserve the subleading single 
In A'' term at O(a^) in Sstan- Indeed, Acoo is already fixed (given coo) from the input 
^00^ = 7oo"' (compare e.g. eq.( |5.3D and ( 5.12| )). Then, the knowledge of Acn fixes A^i 



requiring 

-Po Acii + A^i coi = 0, (5.18) 
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which fixes Ai*^^, given Ai- 

Furthermore one gets cio cio + Acio = Cig"' with 



Acio= / — AG(e)lne. (5.19) 



new 



Thus Acio is fixed (given cio and A^"") from the input -/3o + ^r"'c[Jo'" = 7io 
(compare again eq.( |5.3| ) and ( ^.121 )). One thus sees that in this process AG(e) is determined 
uniquely in principle by its input logarithmic moments Aqq (i : ^ oo) (see eq.( 5.16| ) and 



( |5.19 )), as can be easily checked in the first few orders of perturbation theory. To determine 



A(jr(e), we need however input all the 7i[f""s and use all orders of perturbation theory. This 
procedure is clearly not tractable, so we next turn to large A'^^ where the problem can be 
solved in closed form. 



6. The asymptotic large N Borel transform 

Before doing this, it is useful to first introduce a simplification appropriate to the large N 
limit, valid also at finite Nf beyond the single dressed gluon approximation. I observe that 
the perturbative series eq.(5^) of the Sudakov exponent eq.(5/7) still contains 0{1/N) 
terms at large A^. It is useful to discard them, and find a Borel representation of the 
corresponding series Sas{Q'^ , N) which contains only InA^ and 0{N^) terms, with all the 



0{1/N) terms in eq. (|5.3|) dropped. The (RS invariant) Borel transform of eq.(5/7) takes 
the factorized form (similarly to eq.(4.3) and (4.4)) 

de 

B[Sstan]{u,N) = B[As]{u)exp{ulnN) _G(e) exp(-nlne). (6.1) 

Jo ^ 

Taking straightforwardly the large A'^ limit one gets a tentative ansatz for the looked for 
Borel transform 



/■°° de 

B[Sas]{u,N) = B[As]{u)ex]i{ulnN) / — G(e) exp(-nlne). 

Jo e 



(6.2) 



Since the large Nj calculation is usually performed in the 'dispersive approach' Q which 
makes use of a 'Minkowskian' coupling, it is convenient to use the equivalent 'Minkowskian' 
representation of B[Sas]{u, N), i.e. to introduce the Minkowskian counterpart of the Su- 
dakov effective coupling, together with a 'Sudakov characteristic function' G{y) related to 
the Sudakov distribution function G(e) by the dispersion relation 

where Q = —ydQ/dy. It follows that^ 

Sm TTU poo ^ 

B[Sas]{u,N)=B[As]{u) exp{ulnN) ^^(y) exp(-nlny). (6.4) 

T^u Jq y 

^ B[As]{u) is the Borel transform of the Minkowskian Sudakov effective coupling. 
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Comparing eg. (p. 21) with eg. ( p. 4]) yields the standard relation between 'Euclidean' and 
'Minkowskian' Laplace transforms 







°° de ^, , / , N sinTTU f°° dy ■ / , \ 

— G(ej exp(— u in e) = / — ylyjexpl—ulny). 

TTu Jq y 



(6.5) 



However, eg.( |6.2[ ) and ( |6.4[ ) are not guite correct. For u = the integrals on the right hand 
side are UV divergent, since ^(oo) = G(oo) = —1 (see eg.( ^.6|) ). One actually needs to 
'renormalize' these 'bare' Borel transforms, i.e. to add a subtraction term. One can show 
that the correct answer, in the case of the Euclidean representation, is 



B[Sas]iu,N) = B[As]{u) 



f°° de 

exp(ulnA^) / — G(e) exp(— li In e) 



G(oo) 



u 



whereas for the Minkowskian representation (as follows from eg.(|6.5|)) 



(6.6) 



B[Sas]{u,N)=B[As]{u^ 



sm vru 



exp(nlnA^) / — cj(y) exp(— ulny) ^ 

y 



sm TTU 



(6.7) 



The corresponding 'renormalon integral' representation of Sas{Q , analogous to eg.(|5.7|) 
is 



foo ^r,2 poo JT.2 

SasiQ^N) = —G{Nk'/Q^)As{k^)-G{^) j^^ j^Asi^) (6.8) 

where the second integral on the right hand side on the first line provides the necessary 
subtraction term (the first integral being UV divergent). 



7. Exponentiation of constant terms at large A^^: the single dressed gluon 
result 

7.1 DIS case 

Specializing now to large Nf, the result for the Borel transform oi dlnF2{Q^ , N)/dln Q'^\sdg 
at finite N can be given in the 'massive gluon' or 'single dressed gluon' (SDG) dispersive 
Minkowskian formalism [pO| , |l| as (the V-scheme is assumed) 



B[dlnF2iQ^N)/dlnQ%DG{u,N)=4CF^^-^^ H —^sdg{x, N) exp{-ulnx), 

(7.1) 

where the 'characteristic function' J^sDci^, has been computed in [||] (x = X'^/Q^ where 
A is the 'gluon mass'). I have checked that here too a similar scaling property holds at 
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large N, namely, putting GsDciu, = ^sdg{x, with y = Nx = NX^/Q"^, one gets for 
— > cxD at fixed y 



GsDciy, N) GsDciy, oo) = Qsnciy), (7.2) 

and one can derive from the result of Q in the DIS case that 

QsDciy) = -1 + exp(-i/) - ^ y exp(-y) - ^ y r(0, y) + ]^y^ r(0, y). (7.3) 
Thus for N ^ oo 



S[(ilnF2(g2,A^)/dlnQ2]- (^,iV) =4Cpi_!l!^exp(ulniV) ^^^^^^(y) exp(-nlny), 

(7.4) 

which has also to be subtracted 



B[d\nF2iQ',N)/d\nQ%%aiu,N) = ACp 



1 sinvru 



exp(nlniV) / — ^si:)G(y) exp(-uln?/) + 
y 



where Tsdg{^) = 1 is finite. One can show that 

Tsdg{u) r dy 



u 



y 



u{y) + 1] exp(-ulny), 



(7.5) 



(7.6) 



where i/(y) is the (universal) virtual contribution Q to the massive gluon characteristic 
function for space-like processes^ 



,(y) = - / dz^ ^In 

Jq z-y 



(7.7) 



It follows that 



'^SDg{u) 



TTU 



(7.8) 



^sinvrti/ (1 — n)(l — ti/2) 

The above mentionned fact that i>{y) is universal suggests that Tsdg{u) might also be 
universal for all space-like processes. 

I note that the constant terms of the SDG result are obtained by setting = 1 



in eq.(7^). Thus, if one wants to select as input an arbitrary subset of constant terms 
to be included into a new asymptotic Sudakov exponent SasiQ'^ , ^), one should sim- 
ply change the subtraction function Tsdg{u) in eq.(^), namely define a new input 
B[dlnF2{Q\N)/dlnQ'']%chY 



The normalization is half that of 
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B[d\nF2{Q\N)/d\nQ'fii,a{u,N) = ACp 



1 sin?™ 



/3o vru 



exp(MlniV) / ^e?sDG(?/)exp(-ulny) + '^^^ 

y 



(n) 



(7.9) 



where rg^(^(n) (still with rg'^(5(0) = 1) takes into account the new selected set of constant 
terms to be included together with the InA^'s into the new Sudakov exponent. The latter 
(output) should thus be given by (see eq.(|6.7|)) 



B[Sas]{u,N)=B[A^s"'"]{^) 



sinvrn 



exp{ulnN) j — Q new (y) exp{-u In y) - ^^^^^^ 
y sm iru 



(7.10) 

where Q{oo) = —1 does not change (it determines the leading logs). Since there are no 
0{1/N) terms, one can now identify 



ACFB[Sas]{u,N) = B[dlnF2{Q\N)/dlnQ^]%a{u,N), 
which yields the two master equations 



(7.11) 



B[A^/n{u) r^^„,^(y)exp(-nlny) = -^ T ^^^^^(y) exp(-nlny), 

Jo y Po Jo y 



(7.12) 



and 



B[Ar^]iu) 



1 fsm7ru\ f T^D'ci'^) 



Po \ TTU 

1 ( sin TTU 



It also follows from eq.(|7.3|) that 



/?0 V 



Q(po) 

-pnew ( \ 
J- SDG\^)- 



(7.13) 



/ —GsDGiy)exp{-ulny) 
Jo y 



1 



+ 



(7.14) 



'2 \l-u ' l-u/2^ 

Eg. ( [7. 12 ) and ( 7.13| ) allow to determine both the Sudakov characteristic function Gnew{y) 
and the Borel transform of the associated Sudakov effective coupling B[A^'^^](u) corre- 
sponding to a given input set of exponentiated constant terms which fix the subtraction 
function rg|^(u). It is interesting that B[A'^^'^](u) is given entirely by the subtraction 
term. Once Gnewiy) is determined, one can derive the corresponding Sudakov distribu- 
tion function Gnew{^) using eq.( |6.3D (G„etu(e) is essentially the time- like discontinuity of 
Qnew{y)), thus proving at large Nf the statement made below eq. (|5.9|) . 

In particular, one can easily obtain the result in the two extreme cases where all 
constant terms are included in the Sudakov exponent (the input is then given by eq. (|7.5| )), 
and that where none are. But these two cases are not interesting for the IR finite coupling 
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approach, since one finds that B[A'^^'^]{u) then contains renormalons. In particular, if one 
tries to include all constant terms one gets using eq.( |7.8D 

= (1-.)(Lv2) - 
Despite the presence of IR renormalon poles, this result might be of interest since this 
coupling is expected to be universal for all space-like processes according to the remark 
made below eq.(^i^). However, the corresponding result for the Sudakov characteristic 
function 

, — ^a«Ky)exp(-nlny) = --— ^— ( 1 - ^ ) (7.16) 
y uT{l + u) \ 4 J 

shows that the latter actually does not exist, the inverse gamma function having no Mellin 
representation! 

On the other hand, if no constant terms are included in the Sudakov exponent, one 
finds that the corresponding Sudakov distribution function Giogs{e) is cut-off in the infrared 

Giogs{e) = -eie-l), (7.17) 

so that IR renormalons can show up only through the divergent expansion of the Sudakov 
effective coupling itself. Indeed one gets 

which has IR renormalon poles at u = 1, 2. 

More interesting is the simplest possible ansatz, namely r^*]^^'^(ii) = 1, which gives 
B[A'J""'^']{u) = ^™ and corresponds to ^^'"^'"(A:^) (eq.(|p)), which is already IR 
finite at the perturbative level. However, the first two zeroes at u = 1 and u = 2 lead to two 
log-enhanced power corrections at large in the IR finite coupling framework: the right 
hand side of eq.( |7.14| ) has double IR renormalon poles at u = 1,2, resulting in O(yln^y) 
and 0{y^ In^ y) terms at small y in Gsimpieiu)- In fact the Sudakov distribution function is 
now given by 

Gsimple{e) = GsDcie), (7-19) 

and eq. (|7.3|) indeed yields as e — 



GsDcif^) - - [e{ln€ + -fE - 3) + e^(-lne 



(7.20) 



There is thus a discrepancy with the IR renormalon expectation, in agreement with the 
argument of section 4, and the question arises whether such an ansatz would violate the 
OPE at large already at the level of the two leading (N/Q"^ and N'^/Q'^) power correc- 



tions. Given the well-known |21, 22 1 intricacies of the latter, the question is perhaps not 
yet settled. 
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Applying these results to the standard case, where G{e) is known (eq. (|5.6| )) and fixes 
G{y) from eq.( |6.3[) , and using eq.( [7.12| ), one can also rederive the result eq.(3^), where 



the first two zeroes are cancelled. Here the discrepancy with the IR renormalon prediction 
concerns only the higher order power corrections {N'^/Q^, N'^/Q^,...), which are expected 
22, 10| to be absent at large in the exponent. 



7.2 Drell-Yan case: emergence of an IR finite perturbative effective coupling 

The analogue of eq.(pr^ ) for the scaling violation of the short distance Drell-Yan cross- 



section is 



"^^"Tinly''^^ = i^^y^Q") + s^y(Q"^ ^)) + ^(V^^), (7.21) 



with 

-1 ^A^-l 



SnriQ^N) = / dz 2— AsMi'^ " z?Q^]- (7-22) 

Jo i — z 

Sdy{Q^, N) can again be written as a renormalon integral 

Sdy{Q\N) = —Fr>Y{kVQ\N)As,DY{k^), (7.23) 
where ^5,Dy(^^) is given in eq. (|4.5| ), and the Sudakov distribution function 



FDY{k^/Q\ N) = {1- k/Qf-^ - 1, (7.24) 

involves both [jl3| even and odd powers of k at small k. Taking the scaling limit A — > oo 
with e DY = A^^ / Q fixed one thus get the analogue of eq. (|5.6|) 

Gdy{<^dy) = exp(-e£»y) - 1. (7.25) 

Now the work of for eikonal cross sections suggests to exponentiate a new set of constant 
terms (for any Nf) using 

SiriQ',N)=J^ —G]f^{Nk/Q)A^^^]SY{k'), (7.26) 
with (we deal with the log-derivative of the Drell-Yan cros-section) 

GirieDY) = [Ko{2Nk/Q) + \n{Nk/Q) + , (7.27) 



I.e. 



Gnew r ^ \ 
DY y^DY) 



(IKq 

l + x^{x = 2eDY) 



(7.28) 



where Kq{x) is the modified Bessel function of the second kind.^ But in the Drell-Yan case 
the analogue of eq.( |6.2D is 



GJfy (e_Dy) —1 for eoY oo, consistently with the large env limit of eq.(7.25) 
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B[SDY,as]iu,N) = B[As,DY]iu)exp{2ulnN) / 2— Gz)y(e) exp(-2nlne), (7.29) 

Jo ^ 

and at large Nf we get, instead^ of eq.( 7.12| ) 



f°° df 1 sin TT?; /""^ A; •■ 
B[AT,By]{u) / 2-GSf (e)exp(-2tzlne) = / JigsncDYiy) eM-ulny). 

Jo e m Jo y 

(7.30) 

On the other hand from the large Nj calculation of [14| one gets the A^-dependent part of 
the large Borel transform 

sinyrn dy - t ^ \ 1 r(l - n) 

/ — y5DG,Dy(y)exp(-ulny) = — — — -, (7.31) 

TTU Jo y UT[1+U) 

whereas eq.( 7.28 ) yields 

f °° df 

/ 2— G^ff (e)exp(-2nlne) = -n[r(-n)]2. (7.32) 
Jo e 

From eg. ( [7.30 ) one thus derive the large Nf result 



B[AT:i^y]{u) = -^T^TT— . - 7^^, (7.33) 

Po r(i + ujr(i — u) pq TTU 



i.e. Ag^'^yik'^) is nothing but the effective coupling A'^^"^^'''^ (k"^) of eq.(3.7). The latter, as 
we have seen, does have an IR finite fixed point, and no non-perturbative modification is 
a priori necessary in this case! As expected from the general discussion in section 4, the 
zeroes in S[^^'^|^y](n) lead to large A^ logarithmically enhanced power corrections in the 
IR finite coupling framework, at variance with the IR renormalon expectation. Indeed the 
new Sudakov distribution function has logarithmically enhanced contributions for €dy 



GoYi^DY) = e^y (2 In ez)y + 2-fE - 1) + eV(lneDy + lE-^) + 0(eV In ez)y). (7.34) 

The fact that this new Sudakov distribution function implies a new Sudakov effective 
coupling is of course one of the main point of the present paper. 

8. Reconciling the IR renormalon and the IR finite coupling approaches: 
a large Nf ansatz for a 'universal' Sudakov eff"ective coupling 

I next turn to the question raised in the beginning of section 5 how to reconcile the IR 
renormalon and IR finite coupling approaches to power corrections. Actually, I should first 
stress it is not yet clear whether the two approaches should be necessarily reconciled. For 
instance, in the DIS case, it could be that the OPE at large A^ is consistent with the exis- 
tence of two leading log-enhanced power corrections, as predicted by the Ag"^^'''^ {k"^) ansatz 



Recall also eq.(6.5) 



-17- 



^)), or with the existence of higher order power corrections {N^/Q^, N'^/Q^,...) at 
large in the exponent (see section 7). If this turns out to be the case, the IR finite 
couphng approach would be consistent with the OPE (and at odds with the IR renormalon 
prediction) with the 'simple' ansatz of eq.( |3.7[ ), or eventually with the standard result of 
eq.( |3.4D . Similarly, the 'simple' ansatz of eq.( |3.7D might be the correct one in the Drell-Yan 
case, arising from a 'natural' exponentiation of some 0{N^) terms, despite it also contra- 
dicts the standard IR renormalon expectation. It is interesting to note at this point that 
eq.( |7.13| ) indicates that zeroes in B[Ag^'^]{u) arise from two distinct sources: either the 
'universal' simru/TTu factor (simple zeroes at integer u can come only from there), or the 
'arbitrary' ^^^j^qIu) subtraction term (zeroes at half-integer u can come only from there). 
The previous discussion suggests that zeroes coming from the 'universal' simru/iru factor 
need not be necessarily removed in the IR finite coupling approach, at the difference of 
the more 'artificial' zeroes (such as those occuring in the standard Drell-Yan case eq.( |4.6| )) 
coming from the subtraction term. 

Notwithstanding the above remarks, I shall adopt in this section the attitude that the 
IR renormalon and IR finite coupling approaches to power corrections should always be 
made consistent with one another. For this purpose, one must remove all ^ zeroes from 
B[A'^^^]{u). The mathematically simplest ansatz suggested by eq.(7.13) is to choose 



which yields 



nDGi^^) = ni-u) (8.1) 



Po \ TTU J (3o r(l + u) ' 

It is interesting to compare this ansatz with the result obtained by eliminating the half- 
integer zeroes from the Drell-Yan standard result eq.( f4.6D . There the simplest ansatz is to 
define 

B[AT,lSY]iu) = ^^^11~''^ B[As,dy]{u), (8.3) 



which yields again the 'universal' ansatz eq.( |8.2D . 

Given an ansatz for B [A^^""] (n) , eq. (|7.12| ) then determines Gnew (y) , hence the Sudakov 
distribution function Gnew{^)- In particular, in the DIS case one gets assuming the ansatz 
eq.(PD 



Gne^{e) = -^^^^ < e < 1 



-1 e > 1, (8.4) 



9 



Except those eventually also present in ^Qsoaiy) exp{—ulny), such as the one at it = 4/3 (see 



eq.(|7.14D) 

Apart from the obvious choice "f^s'o'Gi''^) — Trw/sinTru, where the Sudakov effective coupling is just the 
one-loop couphng, and is thus of little interest (having a Landau pole) for the IR finite coupling approach. 
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where e = Nk'^/Q'^. Thus all power corrections beyond the two leading ones are indeed 
absent from the Sudakov exponent, in agreement with the renormalon argument. 
Infrared behavior: although I have not been able to find a closed form analytic expression 
for the Sudakov effective coupling corresponding to eq.(|87 



Arik') = jJ^ <i„exp(-„l„^) (8,5) 

there is very strong numerical evidence that it blows up very fast (but remains positive) 
for A;2 ^ 0: 

Ar(*') = ^exp(|!). (8,6) 

Assuming that eq.( |8.6| ) is correct, one can again speculate that there exists a non-perturbative 
modification SAs^Npik"^) of the coupling at very small momenta which will eventually turn 
the too strong IR singularity into a non-perturbative, but finite IR fixed point. A sim- 
ple ansatz for 5As^i^p{k'^), localized in the IR region, can be constructed in analogy with 
eq.( 3.13| ). Expanding the righthand side of eq.(^.6|) for large fc^ 



1 /A2\ 1 / A2 1 A^ 1 A6 



suggests to subtract term by term using 



1 

To 



A2 

fc2 



exp 



'A2 



1 A^ 



A4 



1 A6 



"A6 



Then for fc^ — > one gets 



SAfS.ik-, . _i_ + + i|! + ...) +1 (l + i + i + ,..) +0(A-VA^ (8.S) 
Thus defining 

AT:^p{k') = A^^^ik') + SA^^^pik') (8.10) 
yields indeed for k'^ a finite and positive IR fixed point 

A-s^^p{k')c.^ + 0{k'). (8.11) 

9. Resummation of the 'left-over' constant terms at large Nf, and a strat- 
egy for phenomenology at finite Nf 

For a given choice of the exponentiated constant terms to be included in the Sudakov 
exponent Snew{Q'^-,N) (eq.( |5.9| )), there remains corresponding 'left-over' constant terms 
HnewiQ"^) (see also eq.(|3T|)) 
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= 4<^^ {Hnen,{Q^) + Snevj{Q^ N)) + 0{1/N), (9.1) 



with 



such that 



7io +A7io =51, (9.3) 
where the gfj's are the fuU constant terms at large 



"^^"^11!$'^^ = + 5o)a.(g2) + (^^^^2 ^ ^^^^ + 5i)a2(Q2) ^ _ ] + ^^^/^^^ 

(9.4) 

From eq.(7^) and ( |7.9D it follows that the Borel transform of HnewiQ"^) at large Nf is given 



by (in the V-scheme) 

PO \ ITU J U 

It is natural to try resumming the HnewiQ'^) series with a renormalon integral representa- 
tion 

HneUQ^) = -^Go,„e»(A:VQ')^r"'(^'), (9-6) 
using the same effective coupling A'^^'^^k'^) as in the Sudakov exponent, which defines the 



'left -ov6r coiists-iits' distribution function Gq new 

{k'^/Q'^). Eq.(^^ yields a representation 

for the Borel transform 

/■oo J 

B[HneAQ^)]{u) = B[A^/^]{u) —Go,nev.ix)exp{-ulnx). (9.7) 

JO ^ 



Comparing eq.(^^) with eq.(9^) allows to determine Go,ne«)(^^/Q^)- For instance, in the 
case of the 'new Sudakov effective coupling' eq.( |8.2| ), one finds that (x = k^/Q^) 



Go,new{x) = 2x [(2x - 1) exp(l/x)r(0, 1/x) - 2x{lnx- je) - (Inx - 7^ + 2)] - ^(x - 1), 

(9.8) 

and one can check that Go^newix) = C(lnx/x) for x — > 00, while for x — > 

G'o,ne«;(x) ~ -2x(lnx-7B + 2). (9.9) 

Eq. (^]^) predicts in the IR finite coupling approach a log-enhanced O (In /Q^ ) A^-independent 
power correction arising from the left-over constant terms. This is consistent with the IR 



renormalon expectation as can be checked by setting = 1 in eq.(7.5) 
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Application to phenomenology at finite Nf-. I suggest to use the same large A'^^ Sudakov 
distribution function Gnewi^) and 'left-over constants' distribution function Go^newix) at 
finite Nf, while adjusting the corresponding Sudakov effective coupling A'^^'^{k'^) to the 
finite Nj situation. For instance, assuming the ansatz eq.( |8.5| ) is the correct one at large 
Nf, one would use the finite Nf resummation formula 



^'",^g-^> = 4C. f '-^G,^^iNW)A'rit^) (9.10) 

+ ACf ^Go,ne^{k''/Q^)Af{k'') + ACf AHneUQ^) + 0{l/N), 

where the {N f -independent) distribution functions Gnew{Nk'^ /Q'^) and Go,nei«(^^/Q^) are 



the same as in eq.(8^) and ( |9.8D , while A'g^"'(k'^) can be determined order by order in per- 
turbation theory at finite Nf in a standard way, matching eq.( |9.10| ) with the perturbative 
expansion of the left hand side (only the A-dependent terms, i.e. the first line, are neces- 
sary for this purpose). I note that at finite Nf there may still remain 'relic' constant terms 
contained in the function /S.HmwiQ'^)^ not accounted for by the two renormalon integrals 



HneUQ^)= -^Go,neUkVQMT'"ik^) + ^HneUQ^). (9.11) 



on the right hand side of eq.( |9.10| ), such that (see eq.(9.1)) 

foo ^1.2 
10 

The expansion of AHmwiQ'^) is expected to be better convergent than that of HnewiQ"^), 
and can be dealt with in renormalization scheme invariant way by using the method of 
effective charges [^]. 

The right hand side of eq.(9.10) is entirely perturbative. To regularize the renormalons 



integrals in the infrared and deal with non-perturbative phenomena (such as power correc- 
tions), one can now assume a non-perturbative modification of the coupling at finite Nf, 
analogous to eg. ( ^.10 ). Its form should be determined phenomenologically. For instance, 



in the spirit of the method of effective charges, one could try a two-point Fade interpola- 
tion between the (known in perturbation theory) weak coupling UV form, and the strong 
coupling IR form of the non-perturbative beta function of the effective coupling A^^'^p{k'^), 
assuming a non-perturbative IR fixed point (its value A^^'^pik'^ = 0) is then one of the free 
parameters to be fitted). A^^'^p{k'^) could then be determined at all scales by integrating 
its own Gell-Mann-Low like renormalization group equation, and reported into eq.( 9.10| ). 



This method (to be investigated in the near future) does not require the introduction of 
any IR cut-off, and represents an interesting alternative (especially if universality of the 
Sudakov effective coupling does hold), with a simple physical interpretation, to the shape 
function approach, as already remarked in section 2. 



10. Conclusions 

The IR finite coupling approach has many attractive features in the context of Sudakov 
resummation. On the phenomenological side, we have seen it is able to yield a very simple. 
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yet non-trivial prediction for the end tail of the Sudakov peak, which goes beyond short 
distance physics. The parametrization of the low momentum piece of the Sudakov effective 
coupling represents also an attractive alternative to the shape function approach. On the 
more theoretical side, I have pointed out the rather unusual occurence of a strongly IR 
divergent, but nevertheless causal, effective coupling, in an all order, but still perturba- 
tive, framework. The essential differences in the IR domain between the standard Sudakov 
anomalous dimensions (such as cusp), which exhibit completely unphysical behavior at 
low scales, and their specific combinations eq. (|3.2| ) and ( |4.5| ) called here 'Sudakov effective 
couplings', where these pathologies cancel out, has been stressed. Only the too strong IR 
divergence localized at the origin requires the introduction of a non-perturbative modifica- 
tion. I feel this case is basically different from the more familiar one of unphysical Landau 
singularities at finite distances, and makes more plausible the (non-perturbative) IR finite 
coupling hypothesis. The simplicity of the proposed non-perturbative large Nf ansatzes 
represents, I believe, further encouragement to support this speculation. Moreover it was 
found that the 'natural' resummation of a set of 0{N^) terms within the procedure of Q 
for eikonal cross sections remarkably leads in the Drell-Yan case at large Nf to the simplest 
example (eq.( p.7D ) of an IR finite effective Sudakov coupling within perturbation theory 
itself, with no a-priori need of a non-perturbative modification. 

The IR renormalon and IR finite coupling approaches to power corrections are po- 
tentially in confiict with each other. We have seen it is possible to make consistent these 
two approaches, by appropriate use of the arbitrariness of exponentiated constant terms 
in Sudakov resummation. This freedom in redefining the Sudakov exponent bears some 
connection with alternative forms of Sudakov resummation previously discussed in the lit- 
terature The latter however involve an infrared cut-off^^, so that IR renormalons 

do not appear through integration over arbitrarily small momenta, but rather through 
divergences of the redefined anomalous dimensions (or, as in [^], the divergence of some 
initial condition sitting outside the exponent). This is in sharp contrast with the present 
proposal, where integration at arbitrarily small momenta is essential, and makes sense 
through the notion of IR finite effective coupling. The mathematically simplest solution 
to resolve the above mentionned conflict leads to the proposal at large Nf of a 'universal' 
non-perturbative ansatz for the Sudakov effective coupling (eq.( 8.10| )). 

Alternatively, it could be that the correct resolution of the confiict favors one approach 
over the other, and that the answer provided by e.g. the 'simple' IR finite ansatz eq.(|3.7|) 
turns out to be the correct one (as suggested by its natural occurence in the Drell-Yan 
process), despite being at odds with the IR renormalon prediction. These issues should 
be resolved by a better understanding of OPE at large N in the DIS case. Another 
related question to be settled is the application of the method of Q to DIS, which should 
determine the corresponding 'natural' Sudakov distribution function and effective coupling. 
Even if the 'simple', purely perturbative IR finite effective coupling ansatz turns out to 
represent an acceptable answer at the perturbative level at large Nf, there may be extra 
non-perturbative effects, similar to those which must occur in the case of the 'universal' 

^^The 'purely logarithmic' Sudakov exponent which excludes all constant terms gives an example at large 
Nf of an alternative form of Sudakov resummation with a sharp IR cut-off, as we have seen in section 7. 
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coupling eq.( 8.10| ), which could also take the form of a non-perturbative modification of the 



effective coupling at low momenta. Whatever the correct choice of the Sudakov distribution 
function turns out to be, it remains for future phenomenological work to determine the 
corresponding form of the Sudakov effective coupling at finite Nf for each process, and test 
for eventual deviations from universality. 
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